BCS-BEC crossover and the disappearance of FFLO-correlations in a 
spin-imbalanced, one-dimensional Fermi gas 
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We present a numerical study of the one-dimensional BCS-BEC crossover of a spin-imbalanced 
Fermi gas. The crossover is described by the Bose-Fermi resonance model in a real space repre- 
sentation. Our main interest is in the behavior of the pair correlations, which, in the BCS limit, 
are of the Fulde-Ferrell-Larkin-Ovchinnikov type, while in the BEC limit, a superfluid of diatomic 
molecules forms that exhibits quasi-condensation at zero momentum. We use the density matrix 
renormalization group method to compute the phase diagram as a function of the detuning of the 
molecular level and the polarization. As a main result, we show that FFLO-like correlations disap- 
pear well below full polarization close to the resonance. The critical polarization depends on both 
the detuning and the filling. 



I. INTRODUCTION 

Ultracold atoms provide a unique opportunity to study 
basic many-body problems both in equilibrium and in 
non-equilibrium situations [l|. A particularly appealing 
feature of these systems is the possibility to change the 
interaction strength over a wide range via Feshbach reso- 
nances. In a two-component Fermi gas, this allows one to 
study the crossover from BCS-pairing to a Bose-Einstein 
condensate (BEC) of strongly bound molecules [2-11]. In 
a situation, in which the two states involved in the pair- 
ing are equally populated, this is a smooth crossover. By 
contrast, in the case of an imbalanced gas, unconven- 
tional superfluid ground states such as the Fuldc-Fcrrcll 
[i[ or Larkin-Ovchinnikov Q (FFLO) state with finite- 
momentum pairs, a Sarma phase with two Fermi sur- 
faces 0, or a mixture consisting of a BEC of strongly 
bound pairs and a Fermi gas of unpaired atoms have 
been proposed 04111 • Experimentally, spin-imbalanced 
two-c omp onent Fermi gase s have first been realized at 
MIT [TJIll and Rice [H 01]. From the spin-resolved 
density profiles and, in particular, the existence of a lat- 
tice of quantized vortices in a rotating gas |l3| , it is pos- 
sible to observe the disappearance of a conventional su- 
perfluid in the center of the cloud with increasing im- 
balance. Assuming that a local density approximation 
applies, this allows one to determine the breakdown of 
BCS-typc pairing beyond a critical imbalance p 3D that 
is close to p 3D ~ 0.4 for the uniform gas at unitarity in 
three dimensions [13, [H| • 

Unfortunately, in the three dimensional (3D) case and 
in the unitary regime, where the scattering length is much 
larger than the average intcrparticle spacing, it is dif- 
ficult, both experimentally and theoretically, to estab- 



lish unambiguously the existence of phases with uncon- 
ventional pairing that are expected when the balanced 
(p = 0) superfluid becomes unstable. The experimentally 
observed density profiles [TtJ at the unitary point are con- 
sistent with the prediction of a first order transition from 
a balanced superfluid to a normal state, in which the two 
spin components each form a Fermi liquid @ ■ This theo- 
retical prediction is based on a variational ansatz for the 
ground state [H, 0j| , which excludes unconventional su- 
perfluid phases. It is therefore of considerable interest to 
study models, for which the phase diagram of the imbal- 
anced gas along the BCS-BEC crossover is accessible by 
methods that are sensitive to states with complex order. 

In the case of one dimension, such powerful numeri- 
cal and analytical tools are indeed available. In fact, for 
both the attractive fermionic Hubbard model [20[ and 
the associated continuum model [U, [22j , there is an ex- 
act solution that can be extended to the imbalanced case 
[23M27I . The ground state phase diagram consists of three 
phases: a balanced superfluid, a polarized intermediate 
phase and a fully polarized, normal Fermi gas [28| . In the 
weak coupling limit, both a solution of the Bogoliubov 
de Gennes equations [29| and bosonization [3(| indicate 
that the polarized intermediate phase is an FFLO-like 
state at any finite imbalance. This prediction has been 
recently verified by density matrix renormalization group 
(DMRG) and Quantum Monte Carlo (QMC) cal- 

culations 36, 37j. It applies both to the continuum case 



and in the presence of an optical lattice, and the FFLO 
state exists in mass- imbalanced systems as well (38l - l4l1 |. 
Moreover, the one-dimensional (ID) FFLO state is also 
stable in the inhomogeneous case that arises in the pres- 
ence of a trapping potential [3ll . |32| . [37j ■ It is important 
to point out that these methods give access to the regime 
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of strong interactions as well, where the energy scale of 
the superfiuid states is of the same order as the Fermi 
energy. In the context of cold atoms, this is the relevant 
regime because in weak coupling, nontrivial order only 
appears at unobservably low entropies of s ~ T c /Tp <C 1 
per particle. 

As realized by both Fuchs et al. [Hf and Tokatly 
however, attractive fermion models are not sufficient to 
account for the full physics of the BCS-BEC crossover in 
one dimension. Indeed, in the strong coupling limit, they 
describe a Tonks-Girardeau gas of dimers. They are un- 
able, therefore, to cover the regime of weakly interacting 
bosons that is reached when the size of the two-particle 
bound state is smaller than the oscillator length of the 
transverse confinement. In this limit, the hardcore con- 
straint of the tightly bound dimers becomes irrelevant. 
Moreover, in models of attractively interacting fermions 
there is only one phase at a finite spin imbalance below 
saturation, namely the FFLO phase E3]. 
As we shall emphasize in this work, the generic phase 
diagram of a more general two-channel model is much 
richer, in particular, close to resonance. 

A description of the ID BCS-BEC crossover that prop- 
erly accounts for the coexistence of fermions and bound 
pairs in the imbalanced case can be achieved in the frame- 
work of the Bose- Fermi resonance model [3, IH| m which 
two fermions in an open channel couple resonantly to a 
diatomic molecule in a closed channel. The associated 
amplitude due to the off-diagonal coupling between the 
open and closed channel determines the intrinsic width of 
the Feshbach-resonance In a continuum description, 
the ID Bose-Fermi resonance model has been studied by 
Recati et al. [46| for the special case of a vanishing im- 
balance, where a smooth BCS-BEC crossover occurs. Its 
BCS side is described by attractively interacting fermions 
while on the BEC side, one has a repulsive Bose gas of 
dimers. In the limit of a broad Feshbach resonance, the 
transition between the two regimes is sharp, yet continu- 
ous. In particular, the quasi-long range superfiuid order 
of the ground state does not change along the full BCS- 
BEC crossover. As realized recently by Baur et al. [13] 
in a study of the associated three-body problem, how- 
ever, the situation is more complex and interesting in the 
case of an imbalanced gas. There, FFLO-physics with 
spatially modulated pair correlations that arc present on 
the BCS-side of the crossover must disappear at a crit- 
ical point, giving room to a Bose-Fermi mixture that is 
a conventional superfiuid, where quasi-condensation ap- 
pears at zero total momentum. At the three-body level, 
this critical point shows up as a change in the symmetry 
of the ground state wavefunction |47j . 

As for studies on the many-body physics of the ID 
Bose-Fermi resonance model, we refer the reader to 
Refs. [4(1 l48M5l| . Bosonization has been applied to the 
balanced case in Refs. [H, |49|, and Bethe ansatz re- 
sults for the imbalanced case have been presented in 
Refs. [5(1 [5l|. FFLO correlations, however, have not been 
discussed in either of these studies. 



Experimentally, the formation of molecules in Fermi 
gases that are tightly confined in two transverse direc- 
tions has been demonstrated by the ETH group [HJ , us- 
ing a balanced mixture. The binding energy of molecules 
is finite for an arbitrary sign of the 3D scattering length a, 
in contrast to the situation without confinement, where 
the two-particle binding energy vanishes on the BCS side 
of negative a. 

The objective of this work is to study a spin- 
imbalanced Fermi gas described by the Bose-Fermi res- 
onance model Hamiltonian. We use a real-space repre- 
sentation with a finite, incommensurate filling and map 
out the zero temperature phase diagram by computing 
pair correlations as a function of polarization and detun- 
ing. We find that FFLO correlations [3, [f| dominate in 
a wide parameter range, and we clarify how the presence 
of molecules affects the stability of this phase. Qualita- 
tively, the presence of molecules binds a certain fraction 
of minority fermions into molecules, reducing the overall 
number of pairs in the FFLO channel. As a main result, 
we determine the critical polarization in the crossover 
region at which FFLO correlations disappear, and its de- 
pendence on filling and detuning. Beyond this critical 
polarization and below saturation, the system is a su- 
perfiuid of composite bosons in the molecular channel 
immersed into a gas of either fully or partially polarized 
fermions. As a numerical tool, we employ the density 
matrix renormalization group (DMRG) method [53M55| . 

This exposition is organized as follows. First, in Sec.HIl 
we introduce the model Hamiltonian and discuss its limit- 
ing cases. Further, in Sec. Ill B[ we analytically solve the 
two-body problem. In Sec. IIII1 we present our DMRG 
results for the pair correlations, the momentum distribu- 
tion, and the number of molecules as a function of filling, 
polarization and detuning. We close with a summary and 
discussion in Sec. ITVl 



II. THE BOSE-FERMI RESONANCE MODEL 

A. Hamiltonian 

We use a minimal Hamiltonian for the one-dimensional 
(ID) BCS-BEC crossover [H, [47| in a real-space ver- 
sion, incorporating the kinetic energies of fermions and 
molecules, the detuning of the molecular level, as well as 
the coupling between the fermions and molecules: 
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c i a ^ S a fermionic annihilation (creation) operator act- 
ing on site i, while mj creates a composite boson on site 
i. The boson energy is shifted with respect to that of 
single fermions by an effective detuning v + 3i. It is cho- 
sen such that the energy for adding two fermions or one 
boson, each at zero momentum, coincide at resonance 
v = 0. The amplitude for the conversion of two fermions 
into a closed channel molecule and vice versa is given 
by the Feshbach coupling constant g. For a negative de- 
tuning v < of the molecular level, it gives rise to an 
attractive two-particle interaction g 2 jv < between the 
fermions (46|. Near resonance v ~ 0, this dominates 
any direct background interaction ?7bg between the two 
fermionic species, which is therefore neglected from the 
outset. The hopping matrix elements for fermions and 
molecules are denoted by t and t mo \, respectively. We 
further set t mo \ = t/2, which accounts for the mass ratio 
of 2:1 between molecules and fermions. L is the num- 



ber of sites. Further, 



yielding the num- 



ber of fermions of each species as N a = J2i( n i-<?)> with 
Nf = iV-f + N± and the pseudo-spin index a =f, 4-- The 
only conserved particle number is N = Nf+2N mo \, where 

iVmoi = E^r'); n ?° l = m l m i- We use n = N / L t0 de- 
note the filling factor and p = (N^ — iVjJ /N as a measure 
of the polarization, which we shall also sometimes refer to 
as imbalance. Note that at maximum one molecule can 
sit on a single site, i.e., the molecules behave as hard-core 
bosons. 



B. Two-body problem and spin gap 

1. Scattering amplitude and bound state energy 

In this section, we calculate the effective interaction 
between two fermions that is mediated by the molecules 
at the two-body level. Following the method outlined in 
[46l |. the bound state energy eb > of two fermions is 
determined by the condition 
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where Do(k,uj) is the bare molecular propagator and 
n(fc, ui) is the self-energy of the closed channel propaga- 
tor (as usual, w and k denote frequency and momentum, 
respectively). 

The resulting equation 



tb-v = g 



dk 



1 



_„. 2tt e& + 4f (1 — cos k) 
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admits a unique, real solution eb > irrespective of the 
sign of the detuning v. Of particular interest is the bind- 
ing energy e* = eb{v = 0) at resonance. Except for the 
scale 2t set by the bandwidth, it only depends on the 
dimcnsionlcss Feshbach coupling constant g' = g/(2t). 
For small coupling strengths g' <C 1, it is given by 
e*/(2t) = .g' 4 / 3 /2 2/3 , while e*/(2t) = g' for g> > 1. The 




detuning V' 

FIG. 1: (Color online) Dimensionless binding energy Q vs. 
detuning as computed from Eq. ©. The thick line is Q, — v 1 , 
the asymptotic behavior in the BEC regime (see the text). 
For comparison, the figure includes DMRG results (squares) 
for the spin gap A at g — t and a small density of n — 0.1, 
extrapolated in system size to the thermodynamic limit L — > 
oo. Inset: characteristic energy e* = e(y' = 0) vs g' . 



ratio e*/ (2t) = 1/ (r*) 2 is essentially the size of the bound 
state (in units of the lattice spacing) at resonance. In 
terms of this characteristic length, the condition for a 
broad Feshbach resonance is simply nr* <C 1 [46| . Tak- 
ing e*(g') as a characteristic energy scale, the equation 
for the dimensionless binding energy = e^/e* for an 
arbitrary value of the dimensionless detuning v' = v/e* 
can be written in the form 



/ 4 + e*/(2t) 
fi[4 + Oe*/(2t)] 



CI, 



(3) 



which is easily solvable for the bound state energy f2(z/) 
as a function of the detuning. The definition of guar- 
antees that CI = 1 at resonance, irrespective of the value 
of the Feshbach coupling g'. In Figure [TJ we show the 
dependence of the binding energy Cl(y') on the detuning 
for three values of g/t = 0.1, 0.5, 1. As suggested by the 
preceding discussion, the CI = f2(z/)-curve is practically 
independent of g' . 

On the BCS side, where v' <C —1, one obtains a very 
small binding energy VU = ^4 + e*/(2i)/(2|i/|) < 1, 
approaching y/Ci = 1/|^'| for small values g' -C 1 of the 
Feshbach coupling. In the BEC regime of strongly posi- 
tive detuning u' S> 1, the binding energy 

Cl = v' + (g/e*) 2 /v' + ... 

follows the detuning, i.e.. the energy of the molecular 
state to leading order. As a result, the closed channel 
fraction 



Z = 



deb 
dv 



= 1 - 



(9/e* 



(4) 



is close to one, as expected in the BEC limit. The di- 
mensionless binding energy CI = (r* /rb) 2 determines the 
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size Tb of the bound state normalized to its value at res- 
onance. For n 3> 1, therefore, this size is much smaller 
than the lattice spacing unless g' >■ 1. 

2. Spin gap 

In the previous section, we argued that the binding en- 
ergy n and in particular, e*, are important quantities to 
characterize the ID BCS-BEC crossover on the two body 
level. We next discuss the relation of f2 to the spin gap A, 
which we calculate with DMRG as a function of filling, 
detuning, and the Feshbach coupling. The connection 
between the binding fl and the spin gap has previously 
been pointed out by Orso (23j . 

The spin gap is computed from 

A(L)=E Q (S Z = l)-E a (S z =0), (5) 

where Eq(S z ) is the ground-state energy of a system of 
length L in the subspace with S z = (JV t - N x )/2. We 
then extrapolate the finite-size data for A(L) in system 
size to the thermodynamic limit L — > oo. 

Figured] includes the DMRG data for the spin gap at a 
filling of n = 0.1 and for g = t (squares). Evidently, the 
spin gap coincides with the two-fcrmion binding energy f2 
not only on the BEC-side v' > 1 where this is expected, 
but also far into the BCS regime. Of course, for very weak 
coupling, this agreement must eventually be violated be- 
cause the spin gap A ~ exp [— 7r/ (2|-y|)] depends on the 
filling n. In particular, it is exponentially small in the di- 
mcnsionless coupling constant \~f\ = l/(2n|ajj) <C 1 (ai is 
the effective scattering length in one dimension, see |42|). 
while the two-particle binding energy e& = e*/v' 2 is in- 
dependent of n and vanishes algebraically with the de- 
tuning in this regime. Near resonance, the spin gap is 
identical with the two-particle binding energy in the low- 
density limit nr* -C 1, as shown by Fuchs et al. j4^ |. 
With increasing values of the filling, however, the spin 
gap increases, as is evident from Fig. [5] The many-body 
spin-gap is therefore clearly distinct from the two-particle 
binding energy 

To illustrate this behavior, we display A as a function 
of filling at g = t in Fig. [Ha) and as a function of g at 
n = 0.6 in Fig.^b), both at resonance v = 0. A = A(g) 
at n = 0.6 also grows with the Feshbach coupling g. 

C. Limiting cases of the Bose-Fermi resonance 
model 

To guide the interpretation of our numerical results to 
be presented in the following sections, we find it useful 
to start with a qualitative discussion of the limiting cases 
of the Hamiltonian Eq. (TTJ) in terms of the dimensionless 
detuning v' = v/e* [see also Ref.l47l which uses the more 
standard opposite sign convention for the detuning] . 

(i) The BEC limit, v' ^> 1 - In this limit, all particles 
are bound in the molecular state, i.e., N mo \ = N/2. At 





FIG. 2: (Color online) Spin gap A at resonance v — vs. (a) 
filling n at g = t and (b) vs. Feshbach coupling g at n = 0.6. 
All results are obtained by extrapolations in system size to 
L — > oo. 

filling N mo i/L < 1, this realizes a superfluid lattice gas of 
hardcore bosons, i.e., effectively a Tonks- Girardeau gas 
of molecules. 

Its ground state is characterized by quasi-long range 
order in the one-particle density matrix 

fflf = (mimj) (6) 

in the molecular channel of the form |/0™ ol | ~ x~ x / 2 
(x = \i — j\) [5^]. As the detuning is decreased and res- 
onance is approached, the molecules start to make vir- 
tual fluctuations into fermions. The presence of excess 
fcrmions suppresses these fluctuations, giving rise to a 
repulsion between fermions and molecules which is pro- 
portional to g 1 jv [I?}. Within a continuum model, this 
effective atom-molecule interaction on the BEC side of 
the resonance has been calculated exactly at the three- 
body level by Mora et al. [oT} . They find that the in- 
teraction is repulsive in the regime where the two-body 
binding energy ej, is larger by a factor 2.2 than its value e* 
at resonance. For smaller binding energies, on the BCS 
side, the effective atom-molecule interaction becomes at- 
tractive and also nonlocal, indicating that the picture 
of bosons that can coexist with unpaired fermions is no 
longer applicable |47ll57j. 

It is instructive to compare the regime v' ^> 1 of the 
lattice model studied here to the corresponding contin- 
uum model studied in Ref. (4(| . In the latter case, the rel- 
evant dimensionless interaction parameter 75 = gs/ns 
(riB denotes the density of molecules) can be tuned to 
values small compared to one even in the deep molecular 
limit because gs ~ | et, | 5 / 2 vanishes as the two-particle 
binding energy |e&| becomes very large. As a result, the 
effective Luttinger exponent K("/b) is then much larger 
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than one and one obtains a weakly interacting gas of 
molecules, whose one-particle density matrix p™ 01 decays 

as \p™j° l \ oc x~ {1/2K) with an exponent 1/(2K) that is 
close to zero. In the continuum and for v' ^> 1, there- 
fore, the weakly interacting molecule gas exhibits almost 
true long range order. This regime, however, is not reach- 
able in the framework of the model Eq. (fTJ) , because even 
in the deep molecular limit v' ^> 1 , where the size of the 
two-particle bound state r\> (in units of the lattice spac- 
ing, see the definition of given above) is much smaller 
than one, we still keep only the eigenvalues and 1 for 
the local molecule occupation number n™° 1 = mjmj. In 
reality, however, more than one closed-channel molecule 
could sit on a lattice site in this limit because the lat- 
tice spacing is much larger than r&. Wc shall not further 
discuss or pursue this question in the present work. Con- 
sequently, while we will be able to see the suppression of 
FFLO physics due to molecule formation, which is the 
main focus of our present work, Eq. ([T]) does not de- 
scribe the full BCS-BEC crossover at a finite imbalance 
that should feature a weakly interacting BEC in the limit 
v 1 > 1. 

(ii) The BCS limit, v' < -1 - Here, iV mol ~ 0. Vir- 
tual transitions into the molecular state give rise to a 
weak attractive on-site interaction U — g 2 /v between 
fcrmions. At a finite polarization p > 0, wc thus expect 
FFLO-likc correlations with real-space oscillations in the 
modulus of the pair-pair correlations 



pair ; f t 

Pij — \ c i,T c i4 c J\T c 
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For small polarizations, these correlations are described 
by the sine-Gordon theory whose ground state is an array 
of domain walls, where the superfluid order parameter 
changes by n HH, H(|. For larger polarizations, the 
domain walls merge and the order parameter acquires a 
purely sinusoidal form with a power law decay 



\Pij I 



oc | cos(Qx)|/x 



a(p) 



(8) 



as a function of the separation \i— j\ = x. The associated 
wave vector 



(9) 



is fixed by the density imbalance via the difference of 
the Fermi- wave vectors k-p >a = nNu/L of the major- 
ity(minority) spins. More precisely, as shown by Sachdev 
and Yang |58[ from a generalized Luttinger theorem for 
Hamiltonians of the form ([T]), the difference fcp^ — kp^ of 
the Fermi wave vectors of the interacting system is quite 
generally fixed by the imbalance p as in Eq. ©. While 
the N a are not conserved separately in the case where 
the bosons are condensed, this theorem implies that the 
wave vector of superfluid order in the fermions is given by 
Eq. ©, independently of the detuning, i.e., the strength 
of the interaction. In the notation of Ref. [9[ , the associ- 
ated FFLO state is thus commensurate. 

The exponent a{p) of the power-law decay has a quite 
interesting dependence on polarization and interaction 
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FIG. 3: (Color online) Number of molecules iV mo i as a func- 
tion of the detuning v for: (a),(b) a resonance with g = Q.lt 
and (c),(d) a resonance with g — t. (a) and (c): balanced mix- 
ture p = 0, different fillings n = 0.1, 0.2, 0.4, 0.6, 0.8 (L = 40). 
(b) and (d): Results for different polarizations at a filling of 
n = 0.6 (L = 40). 



strength, first discussed by Yang [30(. At vanishing po- 
larization p = 0, it is fixed by the Luttinger parameter 
K c > 1 of the attractive ID Fermi gas in the charge 
sector via a(p = 0) = 1/K C . In the limit of small po- 
larizations, bosonization gives a(p > 0) = 1/K C + 1/2 
[30| . i.e., a discontinuous jump of a(p) at p = + . This 
dependence has recently been verified in Ref. [34j |. using 
the attractive ID Hubbard model. 



III. DMRG RESULTS FOR THE IMBALANCED 
CASE 



In this section, we present our DMRG results for the 
number of molecules, the pair correlations, the momen- 
tum distribution function (MDF) of both fcrmionic com- 
ponents, as well as the MDF of the molecules, all as a 
function of polarization, and detuning. As a main result 
we show that, while FFLO correlations are present in 
the BCS limit, as the number of molecules increases, the 
FFLO correlations disappear well below full polarization. 
Upon increasing the polarization at a fixed detuning and 
in the crossover regime, the system thus first has FFLO- 
like correlations, and then undergoes two phase transi- 
tions at polarizations p\ and pi. For p\ < p < P2, pair- 
ing at zero momentum coexists with FFLO correlations, 
while for p2 < p < 1, the system behaves as a Bosc- Fermi 
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mixture with only one fermionic component, the major- 
ity spins. Therefore, the large-p phase is divided into a 
supcrfluid of molecules immersed into either a gas of par- 
tially polarized fermions or fully polarized fermions below 
saturation. We further establish that the molecular and 
pair correlations are identical for p < pi in the sense that 
first, they feature instabilities at the same wave vector 
and second, their highest occupied natural orbitals are 
identical. Our results are summarized in phase diagrams 
for g = t/2 and g = t that are presented and discussed 
in Sec. HITCl 



A. Number of molecules 

To identify the crossover region characterized by a fi- 
nite density of both fermions Nf / L > and molecules 
N mo \/L > 0, we first calculate -/V mo i as a function of the 
detuning v at both g = O.li and g = t. The results are 
depicted in Fig. [3j both for p = and several values of 
the filling n [panels (a) and (c)] and p > at fixed filling 
?i = 0.6 [panels (b) and (d)]. 

We see that in the balanced case, the crossover region 
is between — 3t < v < t for g = O.lt and in the range 
— 4t ^ v < At for g — t. Moreover, the increase of A mo i as 
v is moved from the BCS to the BEC side occurs over an 
increasingly wide range of dctunings with increasing den- 
sity n. This is consistent with the result that an abrupt 
change from a purely fermionic system (iV mo i ps 0) to a 
purely molecular one (Nf ss 0) only exists in the low- 
density limit of a broad Feshbach resonance nr* <C 1, as 
discussed previously in Refs. [H, HD). An obvious, but 
important consequence of the off-diagonal Feshbach cou- 
pling g is that the filling nf = Nf/L in the fermionic 
channel depends on the detuning and the Feshbach cou- 
pling, ranging from nf = n in the v' <C — 1 limit to 
?i/ = in the BEC limit v' S> 1. Therefore, the Fermi 
wave vectors fc^-f/j. vary, too. This is consistent with our 
numerical observation from Fig. [2] Sec. Ill Bl that the spin 
gap is a function of v, n, and g. 

The effect of the imbalance at some generic density n 
[n = 0.6 in Figs. [3Jb) and (d)] is to make the window in 
which molecules and both fermionic species coexist with 
comparable densities narrower. In the g = O.lt case, the 
detuning, at which 2N mo \ ps A, is shifted towards the 
BCS regime v < as the polarization increases. 

Figure Eta) shows the number of molecules 2N mo \/N 
as a function of polarization and for several values of the 
detuning v at g = t and n = 0.6. As soon as the line 
N\. = is reached at some polarization p2, no pairing 
of fermions is possible anymore, and we are left with a 
BEC of molecules immersed into a fully polarized gas of 
fermions. This sets an upper limit, well below saturation 
A = Aj-, for the emergence of FFLO-like correlations. 
In fact, in Sec. IIII C[ we shall see that the FFLO regime 
actually disappears well below p2- 

It is further instructive to compare the polarization de- 
pendence of all particle densities, i.e., majority fermions 



A-j-/ A, minority fermions N^/N, and molecules N mo i/N, 
in the crossover region and before resonance v = —t, 
shown in Fig. Etb). The large-polarization region, in 
which N^/N ps 0, is consequently characterized by a lin- 
ear dependence of A mo i and Af on the polarization, with 
the slope being independent of the detuning v. Note 
that from comparing L = 40 and L = 120 sites data, 
we conclude that finite-size effects are negligible for the 
parameters considered. 

To determine p2, we compute the polarization curves 
p = p(h) for a given detuning and filling n, where h de- 
notes an effective 'magnetic field', coupled to the Hamil- 
tonian through a Zeeman-like term 

Hficid = -h(N t - N±) 

that favors a finite imbalance p > 0. 

The results for g = t and n = 0.6 are displayed for 
vjt = -3,-1,0,1 in Fig. E3 For v = -3t, the p(h)- 
curve has no features, and indicates the presence of a 
very small spin gap. At small polarization, p = p(h) 
increases linearly with h, consistent with recent studies 
of the magnetization process of attractively interacting 
fermions [59l |60| . At v = —t, we first identify the pres- 
ence of a large spin gap (identified by 2h c ), and two kink- 
like features at finite polarizations p\ and P2 ■ Essentially, 
at p > 0, the system is a multi-component Luttinger liq- 
uid, and the presence of kinks indicates the disappear- 
ance or appearance of one component. It is thus easy to 
guess that the kink at larger polarizations, i.e., P2 is as- 
sociated with the depletion of the minority fermions, i.e., 
Aj, ps for p > p2- This is consistent with our results 
for the particle densities shown in Fig. 0Jb) and will be 
further corroborated by the discussion of the momentum 
distribution functions (see Sec. IIII B~Tj) . In view of the re- 
sults for the BCS-BEC crossover of the imbalanced Fermi 
gas m 3D (sec, e.g., Refs. [IS El), one might speculate 
about the possibility that phase separation could appear 
also in one dimension. However, we stress that the criti- 
cal fields h\ and \i2 corresponding to p\ and P2 are well 
separated. In particular, a finite-size scaling analysis of 
the fields hi and /12 for v = — t shows that /12 — hi > 
remains finite in the limit of L — > 00. This rules out the 
possibility of a jump in p(h) and thus of phase separation 
in a uniform system. 

The nature of the first kinkpi in Fig.[5jb) will become 
obvious from the analysis of the pair correlations to be 
discussed in Sec. IIII B] As we shall see, below pi, we have 
pairs at a finite momentum (i.e., the ID FFLO state), 
molecules and the two fermionic components, while at 
p > Pi, additional pairs at zero momentum are formed. 
On resonance, i.e., at v = 0, we still identify a kink at 
P2, while on the BEC side (y = t), the polarization curve 
is smooth, with p(h) oc \/h — h c , where the critical field 
h c for the onset of a finite polarization p ^ is in fact 
connected to the spin gap by the simple relation 2h c = A 

This behavior is characteristic for a band-filling tran- 
sition of a single component, which in this case are the 
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FIG. 4: (Color online) (a) Number of molecules 7V mo i as a 
function of polarization at n = 0.6, g = t, for several values 
of the detuning v. The dashed line shows the maximum pos- 
sible AT mo i at a given p at which N± = 0. The polarization 
P2 at which N mo \(jp, v) meets that line sets an upper limit for 
the emergence of FFLO correlations, (b) Density of majority 
N-f/N (solid lines), minority N±/N (dot-dashed lines), and 
molecules 2iV mo i/iV (dashed lines) as a function of polariza- 
tion for v = -t(n= 0.6, L = 120, g = t). 



majority spins. Note that the same square-root depen- 
dence in magnetization curves has been found for a ID 
Bose-Fermi mixture 1511 . 



B. Pair correlations and superfluidity of molecules 

1. Momentum distribution functions for pairs, molecules, 
and fermions 

To address the key questions of (i) the existence of 
FFLO-likc correlations and (ii) their stability against the 
presence of molecules, we compute the momentum dis- 
tribution function of first, pairs (^? a ) and second, the 
momentum distribution function of the molecules (n™° 1 ) 
by taking a Fourier transformation of the real-space data 
for Eq. ([7]) and of the one-particle density matrix of the 
molecules, pf!j° l [compare Eq. ©], respectively. In the 
following we focus on g = t, unless otherwise stated. 

The results for n£ air and n™°' and a filling of n = 0.6 
are shown in Fig. [5] and Fig. [7J respectively. We choose 
three values of the detuning: v = —3t [panels (a)], which 
is on the BCS side, v = — t [panels (b)] in the crossover 
region, and finally v = [panels (c)] on resonance. It 
is instructive to contrast the behavior of these quantities 
with that of the momentum distribution functions of ma- 
jority and minority spins, i.e., rv k , displayed in Fig. [3] 
n£ is the Fourier transform of the one-particle density 
matrix pfj = (c| )(T c j)(T ). 




FIG. 5: (Color online) Polarization vs. field h (n = 0.6 and 
g = t) for: (a) v — —3t, (b) v = —t, and (c) v = 0, t (thick 
and thin solid line, respectively). The dash-dotted line in (a) 
is a fit to p{h) = a(h — h c ) close to h c , while the dotted line 
in (c) is a fit of the numerical data to p(h) — b\/h — h c . h c is 
the critical field for the breakdown of the BEC at p — with 
N « 27V mol . 

Starting with the Fourier transform of pair correla- 
tions, we note that in the BCS limit and as the polar- 
ization is increased, we observe quasi-coherence peaks at 
a finite momentum Q > [see Fig. Uta)]. Yet, these 
peaks are weak and the pairs' MDF resemble the one 
of a weakly interacting two-component Fermi gas de- 
scribed by the attractive Hubbard model [note that the 
finite-Q peak is more pronounced in the molecules' MDF, 
Fig. UJa)]. The rather weak peaks are probably a conse- 
quence of the fact that the pair correlations differ from 
a pure cosine [as suggested by Eq. ([8])]. This is certainly 
the case at small values p <C 1 of the polarization (see, 
e.g., Ref. [29| and the discussion in Sec. Ill C| . 

The position Q of the maximum in n^ air follows kp^ — 
kp^, as we illustrate in the insets of panels (a) and (b) 
in Fig. [6] This, as usual, is a defining feature of the ID 
FFLO state. 

The quasi-coherence peaks are way more pronounced 
in the crossover region, i.e., v = —t, which is of pri- 
mary interest in this work [see Fig- E](b)] . We observe the 
breakdown of FFLO-likc correlations at a finite polariza- 
tion < pl u < 1. This critical polarization p]P is smaller 
than the upper limit P2 discussed above. An emergent 
feature of the pairs' MDF in the crossover region v ~ — t 
is the coexistence of peaks at both Q = and Q > 
at intermediate polarization [see, e.g., the dotted line in 
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FIG. 6: (Color online) Fourier transform of pair correlations 
at g = t and n = 0.6, as a function of polarization, (a) 
v = — 3t, BCS regime; (b) v = —t, crossover region; (c) v = 0, 
on resonance. The insets in (a), (b) show the position Q of 
the maximum in n^ alr vs. polarization p (squares). The solid 



lines in these insets are kFf 
vs. polarization. 



kF-i = irnp. Inset in (c): 



FIG. 7: (Color online) Momentum distribution of the 
molecules at g = t and n = 0.6 as a function of polariza- 
tion, (a) v — —3t, BCS regime; (b) v = — t, crossover region; 
(c) v = 0, on resonance. The insets in (a) and (b) show the 
position Q of the maximum in n™° vs. polarization p. In 
the FFLO region, the MDF of the molecules shows a peak at 
Q = fci?,-]- — &f,4.. Inset in (c): n^o/^sSoCP = 0) (circles) and 
N mo i /Nmoi (p = 0) vs. polarization (L = 60). 



Fig. Eth)]. By determining the polarization at which we 
see pairing at both Q = and Q > [the dotted line in 
Fig. IDJb)], wc find that this coincides with the first kink 
seen at pi in the polarization vs. magnetic field curves 
shown in Fig. [SJb). Therefore, we conclude that the first 
phase transition and thus the boundary of the ID FFLO 
phase in the crossover regime and at p > is the one at 
p = Pi where pairing at Q = starts to contribute, ef- 
fectively adding an additional quasi long-range order pa- 
rameter to the system. We can further define a crossover 
polarization p* > p\ , beyond which the dominant insta- 
bility is at Q = 0. In the example of v = —t shown in 
Fig. [Sib), p* = 1/2. Note that slightly above p* , some 
modulation in the pairs' MDF survives, which shows up 
as a smaller maximum in n£ air at a finite momentum. Fi- 



nally, we note that the FFLO correlations are typically 
enhanced at low densities (e.g., at n — 0.2; results not 
shown here). To summarize, we identify p]P with the 
upper boundary of the FFLO phase, i.e., p]P = p\. 

Right at resonance (y — 0), no signatures of FFLO 
correlations are visible any more, and the momentum dis- 
tribution functions of both the pairs and the molecules 
feature a maximum at zero momentum [see Figs. EJc) 
and |7(c)] . We observe the same behavior on the BEC 
side, v > 0. For illustration, the k = 0-weight in the pair 
and molecular MDFs are shown as a function of polariza- 
tion in the insets of Figs, injc) and [3c). Quite notably, 
n ™=o exhibits features that can be related to the phase 
transitions the system undergoes as p increases. First, 
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the weight discontinuously drops from its p = value, as 
the critical field for breaking up molecules is overcome 
at p = + . Second, rif^ takes a maximum at p2, where 
the system enters into the Bose-Fermi mixture phase at 
p > P2- A similar, yet less significant behavior can be 
seen in the number of molecules, N mo i(p) / N mo i(p = 0), 
which we have included in the inset of Fig. [TJc) for com- 
parison (solid line) [see also Fig. BJb)]. 

An important point that should be emphasized in this 
context is the fact that the respective quasi-condensates 
of molecules and fermions are locked into each other. In- 
deed, they qualitatively show the same behavior concern- 
ing the position of their maxima, as is evident from com- 
paring Figs. [6] and [7] 

We next discuss the MDF of the two fermion compo- 
nents, shown in Fig. [S] In the BCS limit, the MDFs 
feature a sharp edge, reminiscent of a weakly interact- 
ing lattice gas and consistent with the features observed 
in Fig. Eta). As v moves the system into the BEC 
regime, the p = MDFs become quite broad, as ex- 
pected for a strongly interacting system and and for the 
standard BCS-BEC crossover (see, e.g., Rcfs. 0, H|). 
Upon polarizing the system, n\ develops a sharper edge 
[see Figs. EJa),(b), and (c), left panels], as eventually, 
only the majority fermions remain. This is particularly 
evident in the case of v = — t shown in Fig. |U(b): for 
p > 1/2, N± = J2k n t ~ °- Simultaneously, for p > 1/2, 
n| changes from a smooth function seen at p < 1/2 to a 
steep one, since for p > 1/2, there is a single fermionic 
component left. Thus the depletion of minority fermions 
characterizes the transition to the Bose-Fermi mixture 
phase at p > p2- 

2. Natural orbitals 

To render the analysis of the locking effect pril . |48| . |49| . 
[62[ between p^j 11 and pfj° l more quantitative, we com- 
pute the eigenvalues and eigenvectors of the associated 
one-particle density matrices, p^ air and p™ 1 (the eigen- 
vectors are sometimes called 'natural orbitals'). In par- 
ticular, the orbital <fo that is connected with the largest 
eigenvalue according to the Penrose-Onsagcr decompo- 
sition [63j of the density matrix reveals the real-space 
structure of the quasi-condensates [64[. In the presence 
of FFLO-type order, </>o is therefore a nontrivial function 
even for a homogeneous system. The modulus of this 
quantity, i.e., \<po\ is plotted in Fig. [9ja) for n = 0.2 and 
in Fig. Elb) for n = 0.6; in both cases for p = 0,1/6 
and values of the detuning such that the system is in the 
crossover regime. 

Both at p = and in the FFLO phase, the natural 
orbitals of molecules and pairs are fully identical, as has 
been shown for the limit of vanishing polarization in pre- 
vious studies IHIil- Further, in the ID FFLO phase, 
the spin density 

(S?) = (K t ) - (n u ))/2 



(a) g =t, v=-3t, n=0.6, L=40 




momentum k momentum k 



(b) g=t, v=-t, n=0.6, L=40 




momentum k momentum k 



(c) g=t, v=0, n=0.6, L=40 




momentum k momentum k 



FIG. 8: (Color online) Momentum distribution functions 
at n = 0.6, g = t, as a function of polarization [left panels 
a ="f; right panels: a —\]. (a) v — —3t, BCS regime; (b) 
v = —t, crossover region; (c) v = 0, on resonance. Arrows 
indicate increasing polarization p. In the case of panel (b), 
n l k w for p > 1/2. 



follows the real-space modulation of the natural orbital, 
with excess majority fermions residing in the nodes of 
the quasi-condensate (compare Refs. [2f| Hi] for the case 
of the ID attractive Hubbard model). In contrast to the 
behavior of the spin density, the density of molecules fol- 
lows the modulation of the quasi-condensate. In other 
words, the molecular density has its maxima and min- 
ima at the same positions as the natural orbital. We 
should stress here that the presence of features in the 
densities are due to the open boundary conditions used 
in our simulations. In the limit of L — > oo, the density 
and spin profiles will become flat, while the modulations 
can then be detected in the respective correlation func- 
tions (compare Refs. [^33, 68] for the attractive Hubbard 
model) . In the experimentally relevant situation of har- 
monically trapped particles, however, the density profiles 
themselves should have properties similar to those dis- 
cussed here for finite systems with open boundary con- 
ditions, at least in parts of the particle cloud. 
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FIG. 9: (Color online) Natural orbital \<f>o\ for pairs (lines) 
and molecules (symbols) for (a) n = 0.2, v — 0.005£ and (b) 
n — 0.6, v = —t, both at p — 0, 1/6 (circles and squares, 
respectively). The stars are (Sf) for p — 1/6. Inset in (b): 
density of molecules (n™ 01 ) at p = 0, 1/6. (c) Natural orbitals 
in the intermediate phase pi < p < p2 for v = — t, p = 1/2, 
3 = 4 and i = 60. 
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FIG. 10: (Color online) Decay of pair correlations in real 
space for g = t, n = 0.6, v = — t and p = 0, 1/6, 5/6. Symbols 
denote DMRG results from L — 120 sites while the dashed 
line is a fit of f(x) = a | cos(Q x + <f)\/x aW (Ref. © to the 
numerical data. The fit parameters are a, a and <j>, while Q 
is taken from the Fourier transform of the pair correlations. 
The p — 5/6 curve is off-set by a factor of 0.1 for clarity. 



3. Spatial decay of pair correlations 

To conclude our analysis of the pair correlations, we 
show that the pair correlations at n = 0.6 asymptotically 
decay as |p^ air | oc | cos(Qa;)|/a; a , x = \i — j\, in agree- 
ment with predictions from bosonization for the slowest 



decaying contribution to \p?j \ [30j. To that end, we fit 
f(x) =a\cos(Qx + <j))\/x a 

to our numerical data, measuring j away from the cen- 
ter of the system (i.e., i = L/2). Considering that the 
system sizes are not that large, the agreement between 
the DMRG results and the formula from bosonization is 
remarkable [see Fig. [TU]. In the regime, where FFLO 
correlations have completely disappeared, the pair corre- 
lations decay with a power law, as Fig. [10] suggests for the 
example of p = 5/6. Small oscillations are due to an in- 
homogeneous background density of pairs and molecules 
[compare the inset of Fig. EJb)]. 

Finally, we have also verified that at p = and in the 
BEC limit v' ^> 1, our numerical data are consistent with 
a power-law decay of the one-particle density matrix of 
the molecules 

with an exponent of /3 ~ 1/2. 



Note that in the regime p\ < p < P2, the molecular 
and the pair correlations still exhibit instabilities at the 
same wave vectors [see Fig. HJc)] , even though the natural 
orbitals differ in their amplitude. The locking effect (i.e., 
natural orbitals of pairs and molecules with the same 
amplitude) is re-encountered in the high- field region pi < 
p < 1. There, the molecular \4>q\ is smooth, while the 
corresponding natural orbital for the pairs exhibits small 
oscillations. 



C. Phase diagram 

Our results for the phase diagram of the ID BCS-BEC 
crossover described by Eq. ([1]) are summarized in Fig.lTTl 
for the cases of g = t [panel (a)] and g = t/2 [panel 
(b)]. The main panels contain the data for n = 0.6 and 
we present polarization p vs. dimcnsionless v' detuning 
phase diagrams. 
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FIG. 11: (Color online) Phase diagrams, polarization p vs. 
dimensionless detuning v' for n = 0.6 and (a) g = t and 
(b) g — t/2. The line p = Pi(v') (squares) separates ID 
FFLO from the BEC+PP LL regime. The stars denote p2 (see 
Sec. HITA)l . separating BEC+PP LL from BEC+FP FG. In- 
sets in (a),(b): dependence of p\ on filling: n — 0.6 (squares), 
n = 0.2 (diamonds), n = 0.1 (circles), (c) The same as in (a), 
yet here plotted as a magnetic field h' vs. detuning v 1 phase 
diagram. h c (triangles) is the critical field for the breakdown 
of the balanced gas, while h sa ,t (circles) is the saturation field. 
Lines are guides to the eye. 



We identify three regions at p > 0: (i) the BEC limit, 
v' ^> 1 and P2 < p < 1. Here, molecules are immersed 
into a sea of fully polarized fcrmions. This phase is 
denoted as BEC+FP FG in the figures, where FP FG 
stands for fully polarized Fermi gas. 

(ii) The ID FFLO phase at < p < px. In 
the crossover regime, FFLO is suppressed as p is in- 



creased. We have determined the phase boundary pi 
(open squares) from both the position of the first kink 
in the polarization curves and from the pair correlations. 
In the latter case, at pi, the peak at Q = starts to 
build up in the MDF of the pairs. For instance, the ID 
FFLO phase extends up to v < — 0.3£ at this filling and 
g = t. This is slightly before resonance on the BCS side, 
where, nevertheless, the density of molecules is already 
finite, i.e., N mo i > (compare Fig. [3j) . Lastly, there is a 
region (hi) pi(n, v) < p < P2, beyond which we have a 
Q = superfluid of molecules immersed into a partially 
polarized (PP) fermionic gas. This third phase, denoted 
by BEC+PP LL, is eventually replaced by the BEC+FP 
FG phase at p > P2, where we determine pi from the 
analysis of p = p(h) curves (see Sec. MI AT) . 

Note that the boundary of the ID FFLO phase, p\, 
depends on the filling n. From the insets of Figs. [TTTa) 
and (b), we infer that the larger n, the wider the crossover 
region is, consistent with the discussion of the number of 
molecules (compare Sec. IIII B~T|) . As n — > 0, the critical 
line p = pi becomes quite steep and approaches v rj 
(0.048 ± 0.002)i, or u' w 0.97, for g = t. 

The comparison of the g = t and the g = t/2 phase 
diagram shows that the FFLO phase disappears much 
faster in the case of g = t/2, well before resonance. 
Qualitatively, one can ascribe this to the fact that with 
decreasing values of the Feshbach coupling the number 
of molecules, or more precisely, the closed channel frac- 
tion [compare Eq. @] becomes larger. The presence of 
molecules tends to reduce the number of pairs with FFLO 
correlations. This can be expected to more efficiently 
suppress FFLO physics the smaller g is since the lock- 
ing of molecules and pairs is then also weaker. These 
observations are consistent with our DMRG results for 
the number of molecules and their dependence on polar- 
ization and detuning presented in Fig. [3] In particular, 
the maximum number of molecules is reached at smaller 
values of v the larger the polarization is. 

Figure [TlTc) shows the data of panel (a) in the mag- 
netic field vs. detuning plane, using the dimensionless 
detuning v' and field, h' = h/e*. This yields additional 
information on the saturation field h sa ± and the zero-field 
spin gap A of the standard ID BCS-BEC crossover of the 
balanced system, measured by h c . In comparison with 
Fig. [TJ where we have shown A ~ e* for n = 0.1, we re- 
peat that the spin gap A = 2h c is an increasing function 
of the filling n [compare also Fig. Ufa)]. In the limit of 
v' 1, A = 2h c behaves as A tx v' since there, indepen- 
dently of filling, the ground state of the balanced system 
has N w 2JV mo i and N f w 0. 

In a previous work on the three-body problem in the 
continuum limit, Baur et al. |47j have shown that the 
change in correlations between an oscillating behavior 
on the BCS side due to FFLO physics to a smooth one 
on the BEC is revealed in the symmetry of the three- 
body ground state wavefunction. The numerical value 
of the detuning where this change occurs is v' c ss 0.63 
1471. It is remarkable that a similar critical value for the 
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disappearance of FFLO correlations is also found in our 
many-body calculation of the phase diagram. Indeed, in 
the low-density limit, where a comparison makes sense, 
the boundary of the ID FFLO at small polarizations is 
typically close to resonance, yet on the BEC side of pos- 
itive detuning v > 0. For a quantitative comparison, we 
have determined the critical value v' c (n = 0.1) for the 
loss of FFLO correlations for several values of g from 
data taken with L = 120 sites and polarization p = 1/6, 
the smallest imbalance possible for this system size. The 
resulting values are in the range of 0.55 < v' c < 0.91, 
remarkably close to the value inferred from three-body 
physics in Ref. [47 1. 

In conclusion, it is evident from Fig. [11] that the best 
regime for observing the ID FFLO state is (i) low density 
and (ii) small polarizations. The low density will favor 
a large weight in the quasi-coherence peaks, while the 
polarization needs to be kept smaller than p\. Moreover, 
the ID FFLO phase is more stable at large Feshbach 
couplings g. 



IV. SUMMARY AND DISCUSSION 

In this work, we studied the Bose-Fermi resonance 
model in the imbalanced case as a simple model to de- 
scribe the BCS-BEC crossover of a spin-imbalanced sys- 
tem in one dimension. Our main focus was on the exis- 
tence and stability of the ID FFLO phase. So far, many- 
body calculations of ID FFLO physics were mostly con- 
cerned with models of attractively interacting fcrmions, 
which do not account for the existence of composite 
molecules in the closed channel, typically encountered 
in experiments. Using a numerically exact method, the 
density matrix renormalization group method, we com- 
puted several quantities to characterize the crossover, in- 
cluding the number of molecules, pair correlations, the 
momentum distribution function, as well as polarization 
curves. Most notably, we found that FFLO correlations 
are suppressed in the crossover region due to the pres- 
ence of the diatomic molecules. In particular, the ID 
FFLO phase gives room for a regime of molecules, quasi- 
condensed at zero momentum. The latter is first im- 
mersed into partially polarized fermions, which is then 
replaced by a Bosc-Fcrmi mixture with spinless fcrmions 
below saturation. Thus, the system undergoes two phase 
transitions in the crossover region at critical polarizations 
Pi < Vi < 1 as the polarization increases. 

While our work was concerned with the homogeneous 
system, in experiments, the particles typically experi- 
ence a confining harmonic potential. The shell struc- 
ture for attractively interacting fermion models in ID 
was intensely discussed. The emerging picture for the 
continuum case, based on numerically or analytically ex- 
act approaches (the latter ty pica lly combined with the 
local density approximation) [23l [24 . l2rl [13] is that one 
finds either fully paired wings at small polarization or 
fully polarized wings, while the core is always partially 



polarized. In the case of lattice models, DMRG calcula- 
tions that take the trap into account exactly report fully 
polarized wings with a partially polarized core [3lL 1321 ] 
at intermediate and large polarizations, and the latter 
also remains true in coupled chains at sufficiently large 
polarizations [65| . 

While we expect the behavior of trapped, attractively 
interacting fermions to carry over to the BCS regime 
of the Bose-Fermi resonance model, a finite density of 
molecules may lead to qualitatively different shell struc- 
tures. For instance, the heavier molecules should mostly 
reside in the center of the trap. On the one hand, one may 
expect this to destabilize the FFLO phase in the core, 
while on the other hand, as long as the Feshbach coupling 
g and hence, the locking between pairs and molecules is 
sufficiently strong, the locking could protect the FFLO 
correlations. The clarification of the effect of a harmonic 
trap is left for future research. 

An important question is how the FFLO state can be 
detected in an experiment. Several proposals have been 
put forward, for instance, timc-of-flight measurements 
[66| , the analysis of noise correlations J33, [67} , or features 
in the spin density and correlations [68|. Regarding the 
spin correlations, one expects a peak at nonzero momen- 
tum 2Q ^= in the presence of FFLO order [g^. In fact, 
the spin density follows the modulation of the natural 
orbitals, as has previously been demonstrated for the ID 
attractive Hubbard model [3l|. As we showed here, this 
behavior is also realized in the FFLO phase of the Bosc- 
Fermi-resonance model (compare Fig. Even if the 
FFLO phase was present in a 3D system, the obstacle 
there is that, if at all, the FFLO phase is in the wings of 
a 3D, trapped Fermi gas (see, e.g., Ref. [1^|). This con- 
stitutes another advantage of searching for FFLO physics 
in a ID system: there, the core of a trapped gas will host 
this phase [23|, HtJ , and therefore, the associated mod- 
ulation in the spin density should exist in a large part of 
the cloud, contrary to the 3D case. 
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